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Modular Hamiltonian of Excited States in Conformal Field Theory 
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We present a novel replica trick that computes the relative entropy of two arbitrary states in 
conformal held theory. Our replica trick is based on the analytic continuation of partition functions 
that break the replica symmetry. It provides a method for computing arbitrary matrix elements of 
the modular Hamiltonian corresponding to excited states in terms of correlation functions. We show 
that the quantum Fisher information in vacuum can be expressed in terms of two-point functions on 
the replica geometry. We perform sample calculations in two-dimensional conformal held theories. 


In recent years entanglement theory has found numer¬ 
ous applications in the study of quantum phases of mat¬ 
ter, relativistic field theories and gravity. Most of these 
applications focus on an entanglement measure in bi¬ 
partite pure states known as the entanglement entropy. 
Unfortunately, in relativistic field theories entanglement 
entropy suffers from ultraviolet divergences. In gauge 
theories the definition of entanglement entropy is am¬ 
biguous [1]. In this letter, we present a method to com¬ 
pute in held theory another measure called relative en¬ 
tropy that is provably ultraviolet hnite, universal and free 
of gauge ambiguities mm- 

Relative entropy is a measure of distinguishability be¬ 
tween two states and has nice monotonicity and positiv¬ 
ity properties. It appears naturally in the dehnition of 
entanglement measures for mixed states such as mutual 
information and the relative entropy of entanglement [3] . 
Recently, thinking in terms of relative entropy in quan¬ 
tum held theories coupled to gravity has led to new de¬ 
velopments such as a proof of quantum Bousso bound 
[2, and the identihcation of new gravitational positive 
energy theorems [5]. 

The relative entropy of the density matrix cj) with re¬ 
spect to if is dehned to be 

S{(l)\\ip) =tr{(l) log (j))-tr{^ log Ip). (1) 

Note that relative entropy is ill-dehned when ip is pure. 
The relative entropy of two states can be thought of as 
the expectation value of the difference of the modular 
Hamiltonians of the two states 

s{tP\\ip) = imw-Him) 

= tr{{cP-iP)Hm-^S. ( 2 ) 

Here the positive Hermitian operator H(ip) = — log ip 
is the modular Hamiltonian of ip, and AS' is the dif¬ 
ference of the entanglement entropies of (p and ip. If 
we formally dehne the generalized free energy function 
F-ip{(p) = tr{(pH{ip)) — S{(p), then the relative entropy is 
the free energy difference between the two states 

S{Pi\\iP)=F4cP)-F^{iP). (3) 

The function has all the properties one expects from 
free energy in a thermodynamic theory where ip = 






(o) 


( 6 ) 


FIG. 1: (a) Operator-state correspondence in radial quantiza¬ 
tion of conformal field theories. (6) Reduced density matrix 
corresponds to a path-integral with two operator insertions 
and a cut on the subsystem. 


plays the role of the equilibrium state m |B]. Note that 
F.ijj achieves its minimum on the equilibrium state ip.^ 

In this letter, we construct a class of field theory parti¬ 
tion functions that are proportional to tr{(pip'^~^). Their 
analytic continuation provides the relative entropy and 
the modular Hamiltonian of density matrices in excited 
states \(p) and \ip) reduced to the subsystem. While the 
formalism presented here applies to all quantum held the¬ 
ories we specialize to conformal held theories to have ac¬ 
cess to more computational tools. 

According to the operator-state correspondence in con¬ 
formal held theory (CFT) there is a one-to-one map 
between wave-functionals and operators in the Hilbert 
space. In radial quantization, the wave-functional of an 
excited state \ip) is found by performing a Euclidean 
path-integral with the corresponding operator ft in¬ 
serted. Restricting to subsystem A the state is described 
by a density matrix ip a; see hgure[2 To simplify nota¬ 
tion we suppress the subsystem index A, and use ip to 
refer to the reduced state. 

In principle, one can compute the logarithm of den¬ 
sity matrix directly from the path-intergral by taking 
the logarithm of a path-ordered operator using the so- 
called Magnus expansion, however in practice this is too 
hard. Here, we propose an alternative method to com¬ 
pute matrix elements of the modular Hamiltonian of ex¬ 
cited states from the analytic continuation of correla- 


^ This is a consequence of positivity of relative entropy: > 
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the modular Hamiltonian of excited states 


S{(j)\\ip) = dnlog 


- H{ao)\<j}) = log 




J n—^l 

11 ■ 




n—1 


(5) 


FIG. 2: (a) Entanglement entropy replica trick: the Eu¬ 
clidean path-integral on the n-sheeted manifold corresponding 
to the partition function (6) The .Z„-breaking parti¬ 
tion that appears in our relative entropy replica 

trick. 


tion functions. Our method is a generalization of the 
replica trick in laiH] to the case where one breaks the Z„- 
symmetry among replicas. This enables us to compute 
matrix elements of the modular operator for all states. 
This is in contrast with the old replica method which was 
restricted to states that have local modular Hamiltoni¬ 
ans, the only known example of which is vacuum reduced 
to a half-space or spherical subsystems pnu].' 


RELATIVE ENTROPY AND MODULAR 
HAMILTONIAN 


Consider the Hermitian operator = 

(j)) built out of reduced density matrices (j) 
and il) corresponding to global states \(j)) and |■)/’), respec¬ 
tively. Its trace in conformal field theory corresponds to 
the n-sheeted partition function The idea is 

to take advantage of the analytic properties of correlators 
by using the operator identity 


tr((^logV’) = (4) 

Our partition functions of interest, break the 

Zn replica symmetry present in both the Renyi entropy 
and Renyi relative entropy replica tricks mm- In con¬ 
trast with the symmetric case, our partition functions 
are not monotonic in index n, and have no known oper¬ 
ational interpretations. Nonetheless, under the assump¬ 
tion of analyticity, they provide a computational tool for 
finding relative entropies and the diagonal elements of 


^ In two dimensional CFTs, a finite temperature state on a line 
also has reduced states with local modular Hamiltonian m 


where cto is the reduced density matrix in vacuum. We 
subtract the vacuum modular Hamiltonian so that we 
have ultra-violet finite quantities at any n. The off- 
diagonal elements of the modular Hamiltonian are ob¬ 
tained from its diagonal element in superposition states; 
see supplementary material. 

Each of the terms inside the logarithm above can be 
expressed as a Euclidean path-integral with operator in¬ 
sertions on a replicated or the original geometry [13] . For 
instance, consider the terms Sewing n copies 

of the density matrix cyclically along the boundary of 
their subsystems we obtain n-sheeted replica manifold 
TZn and 2n operator insertions (figure 

m 

( 6 ) 

where Zi and z' are points z and z' on the sheet of 
Tin- It is important to note that plugging Q in ([^ all 
partition function terms cancel and we are left only with 
correlation functions at any n which are free of ultraviolet 
divergences. 

Written explicitly in terms of correlation functions we 
find the main results of this section: 


sm^) = 

9nl0g 




/Til 


dn log 


mz'nMzn)o^r'^) 


{x\H{i;)-H{aoM) = dn 


log 


X-i 


L^+iJ 
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X. 
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where 






= + |cp E^Ao^A~^^) 

+cE^A^t~'^^) + h-c. 

(O) = i'^An)x{zn)0)n„ 

V i‘^A[)Hzi))nAx{z[)x{zi))ni 


( 7 ) 


Here, we have assumed that (tpA) = 0- 
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QUANTUM FISHER INFORMATION 

Our replica trick connects the modular Hamiltonian of 
excited states to analytic continuation of 2n-point corre¬ 
lation functions. Apart from integrable models and large 
central charge theories obtaining analytic expressions for 
2n-point functions is an intractable problem. However, 
as we show in this section a great simplification occurs 
once we focus on near-vacuum states. 

Let us first consider a one-parameter family of states 
Perturbed around \(j)) in perpendicular direc¬ 
tion |A). The reduced density matrix on subsystem A 
expanded in e has the form 



Relative entropy is a smooth non-degenerate function of 
two states. Hence, the relative entropy of two nearby 
states expanded in e vanishes to the first order. The 
coefficient of the second order term, F^{X,Y), is called 
the quantum Fisher information at point cj) in the space 
of density matrices: 

Sicl> + epx\\(t>) = e^F4X,X) + Oie^). 

This function defines a metric on the space of perturba¬ 
tions to state (j) 


2F^{X, Y) = F^{X + Y,X + Y)- F^{X, X) - F^{Y, Y). 


Quantum Fisher information is a local measure of dis- 
tinguishability, and is intimately connected with uncer¬ 
tainty relations [T3] . Consider the relative entropy of two 
nearby states. Our replica trick in ([^ implies 


F^{X,X) = 

- 


m—0 


{^^^)nAXX)n, 


where X and $ denotes the operators that create the 
perturbations corresponding to \X) and \4)), respectively. 
The location of operator insertions are the same as 
For near vacuum states, we replace $ in ([^ with iden¬ 
tity. The quantum Fisher information takes the form of 
an analytic continuation of two-point functions on the 
replica geometry; see hgure 


F,{X,X) = dn 


^ Kx{z"^,z ) 


- n—>-1 




This implies that in arbitrary dimensions the vacuum 
Fisher information of any primary excitation reduced to a 



FIG. 3: The type of two-point functions on the replica man¬ 
ifold whose analytic continuation determine quantum Fisher 
information in vacuum. 


ball or radius R is universal in the sense that it depends 
only on energy and subsystem size. In the remaining 
of this letter, we provide examples of relative entropies, 
modular Hamiltonians and quantum Fisher information 
in two-dimensional CFTs computed using the method 
above. 


EXAMPLES IN TWO DIMENSIONS 


Relative entropy of excited states: Consider a free 
massless boson CFT in two dimensions on a circle of 
radius R and a subsystem at A = (—^/2,//2). We are 
interested in the excited states obtained by the action 
of chiral vertex operators on vacuum at past infinity: 
|q^) = Ia|Q) = where is the boson field. The 

dimension of this operator is {h,h) = (a2/2,0). Here 
X = l/R is the dimensionless parameter. In supplemen¬ 
tary material it is shown that in two dimensions one can 
equally use correlators on cylinder, full complex plane or 


a strip in our formulae in (251 for relative entropy and 


modular Hamiltonian. The conformal factors found from 
the change of coordinates vanish in the limit of n —>■ 1. 
In a free theory with a non-degenerate ground state all 
_>i correlation functions are determined by Wick’s theorem 

(8)m- 


where S refers to correlators on a strip of width 27r, and 
9n = nl^(lx]n ) ' holomorphic excitations [14]''' ^ 
V-a- Therefore, 


S(c 


(oi^)s(V-pVp}^ 


= dn log 


n—1 




= (a — /I)^(I — 7ra;cot(7ra;)). 


( 10 ) 


The analytic continuation used above is justihed in sup¬ 
plementary material. When /3 = 0 this matches the re¬ 
sult previously found using a Z„-symmetric replica trick 
in [12]: S'(a||0) = a'^(l — TTXcot{Trx)). Interestingly, the 
answer in is symmetric in its arguments, 5'(a||/3) = 
^(/Illa). These excited states further have the property 
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FIG. 4: (a) The n-sheeted manifold corresponding to the par¬ 
tition function . (6) The 2n,-breaking partition = 

tr{aAB{cTA ® 



FIG. 5: Inserting the resolution of identity in Z^^ we observe 
that at each K we multiply sphere one-point functions that 
are zero unless 4>x is identity. 


that ^(a) = S(/3) = S'(cro), where (Tq is the vacuum den¬ 
sity matrix. Hence, we find tr{paHj 3 ) = tr[ppHa) for all 
a and /3. 

Modular Hamiltonian of excited states: In the free 
c = 1 CFT Wick contractions imply that a correla¬ 
tor is zero unless all a 7 we have 

{V-aVjO^I^~^^)s = 0. As a result, Xc in (25) is indepen¬ 
dent of c, and we find that the modular operator Hp has 
no off-diagonal terms in the |a) basis. 

The diagonal elements are 


of the vacuum partition functions on manifolds and 
illustrated in figure]^ 

I{A : B) = lim ^ (log Z^^ - log . (13) 

n-J-l n — 1 

The first partition function Z^^ corresponds to Renyi 


entropies of aAB- Therefore, from (12) all we need to 
check is 


dnZr, 


A,B\ 


^^=S{A) + SiB). 


(14) 


{a\H{/3) — H{ao)\a) = fi{(3 — 2a)(l — ttcc cot( 7 ra;)). 

Note that in the limit a = P this reproduces —/3^(1 — 
7 rcot( 7 ra;)) = —S'(/3|tTo) as it should. In the limit /3 = 0 
the difference of modular Hamiltonians is the zero oper¬ 
ator and hence the answer should vanish as it does. 

Quantum Fisher metric around the vacuum: Con¬ 
sider an arbitrary two-dimensional conformal field theory 
on a circle. The vacuum Fisher information is given by 
equation (§. After some algebra we find 

FcriX, X) = dn [2s„(0) -b Sn{x) + S„(-x)]„^;^ 

( ■ 2 / N \ n-1 

) n sin ( 7 r(m + 


For simplicity we expand in small x to find: 




^ Tree ^ 


n sm(7rm/n)-2(''+^i 


2(7ra;)2fo+foV^r(/i-b h-b 1) 

r(/i + h+|) ’ 


n—>-1 
( 11 ) 


where we have used the analytic continuation found in 

m- 

Multiple intervals: The replica trick developed here 
can be applied to subsystems with multiple intervals. As 
an example we focus on mutual information in vacuum 


S{aAB\\crA<Z> ctb) = I{A : B), (12) 


where A and B are non-overlapping intervals. According 
to ([^ the relative entropy is the analytic continuation 


Riemann-Hurwitz formula tells us that Z^^ has genus 
(n—1) and Zf^F jg simply the Riemann sphere. Following 
m we compute the path-integral over these manifolds 
using twist operators in an orbifold theory with replica 
copies of the fields. In particular, up to normalization 
is the correlation function 

('^A)^(n---l) (’^A)^(n.-.271—1) (^il)^(2n—l-.-n) (^s)) 

in a (2n — 1) replica theory. Here ua and 
va are the endpoint of interval A, and going 
around the twist operator cr(i...m) the replica fields 
transform as (X^, Al^ • • • X™, —>■ 

(X2,...X'",X\X™+fo--X2’"-i). Inserting a resolu¬ 
tion of identity on sheet splits the correlator into a 
sum over the product of sphere one-point functions, (see 
figure]^. Sphere one-point function are zero unless 
is the identity operator. In other words, 

S:^F = s^ + S^ (15) 


which is the sum of Renyi entropies of intervals A and 
B\ and hence (14) follows. 


DISCUSSION 

In this letter, we have developed a replica trick that 
takes advantage of breaking the replica symmetry to 
access the modular Hamiltonian of excited states. In 
the absence of the Z^ replica symmetry Renyi's are not 
monotonic in n, hence our method cannot be used to 
obtain lower or upper bounds on relative entropy. The 
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applicability of this method crucially relies on our ability 
to analytically continue correlation functions in n. Ac¬ 
cording to the Carlson theorem [17], in order to find the 
unique analytic continuation of Renyi’s at integer n one 
needs to further fix the behavior at n —>■ ±*oo. We post¬ 
pone a careful study of this asymptotic choice and its 
physical implications to future work. 

The correlation functions needed to compute the mod¬ 
ular operator of an excited state are 2n-point functions. 
There are not many examples of CFTs for which we have 
access to high-point correlators. One class of such CFTs 
are free theories which we briefly discussed. Another class 
are CFTs with large central charge, where one can reduce 
the calculation of n-point functions of heavy operators to 
a classical monodromy problem for differential equations 
that correlation functions satisfy m- 

In holographic theories, the vacuum Fisher information 
in spherical subsystems was recently shown to be dual 
to canonical energy in gravity |2Q|. This confirms the 
universal feature suggested by equation It would 

be interesting to understand the connection between the 
CFT calculation of this quantity and canonical energy in 
the bulk. 
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Superposition states 


In the radial quantization, the wave-function of an ex¬ 
cited state is a path-integral over the disk with unit ra¬ 
dius and an operator inserted at the center of the disk: 

{(j)a\<l>)= [ I?(/)(a;)e-^['^l^(u; = 0), (16) 

J (pa 


where the boundary conditions are (j){uj = 1) = (j)°‘. It 
is immediate that a superposition state of j^) and |x) is 
given by 

C^(0a|V^) +Cx((/>a|x) = [ T><))e"^(‘^^(c^«'(0) -bCxX(O)). 

J (pa 


The complex conjugate state {4>\ is created by the action 
of [4>{z,z)]^ = z~‘^^z~'^^(l){\/z,\/z). A simple way to 
keep track of normalization for superposition states is to 
normalize term by term 


(CclTlICc) 


(x^x) c(^^x) 

($$) (xx) \/($$)(xx) 


We are interested in off-diagonal elements of the modular 
operator which are obtained from its diagonal element in 
superposition states 

{x\M4’) = 2 (^( 1 ) ~ ^(- 1 ) + ~ ’( 1 ’^) 


where (^c|A|^c) = A(c) and j^c) = ;y=j^(l<^) + c|x))- 

Note that \(j)) and |x) are orthogonal: (x®) = 0. 

The expectation value of the modular Hamiltonian in 
superposition state is: 


log(^'(^o)^(io))7?,i - 9„log 


A, 


LA$$(1) + |c|2A^^(I)J 


From 0 we find 
(xjFlW --ffolc/-) = 


log 


A i 


lA+ij 


flog 


fA-t-iJln-)-! 
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where 




= + |c|2 

+cE^^{0^^-^'>) + h.c. 

^ {Hz’u)x{zn)0)n„ 

\/(®(4)^(^i))ki(x(z'i)x(^i))ki 


Conformal map from cylinder to complex plane and 

strip 

There are two choice of coordinates that are more con¬ 
venient for computations in two-dimensional CFTs: the 
complex plane coordinates (w), and the strip coordinates 
(t). Our starting point is the cylinder path-integral with 
operators inserted at z = Eioo and boundary conditions 
(j)a and ipb above and below A. We first apply the map 

sin(7r(z -f l/2)/R) 
sin(7r(z — l/2)/R) 

that sends the cylinder TZi to the complex plane. Under 
this transformation, the subsystem A is mapped to the 
half-line (0, oo), and operators insertions are at 
Then, we apply a second conformal map: lj = e**. The 
complex t-plane is a strip of width 2tt with boundary con¬ 
ditions imposed on t = ±7r lines, and operators inserted 
at t = ±7ra;. Similar maps 


where the subscript S refers to the correlators on a strip, 
and 


JsiU)= ('I 


(n —1) A 


dz 

dl 


(n—1) A 


n—1 / o \ —A / a-\ —A 

az\ oz 


’‘rivai,. 

i=l ^ ^ * 


dt 


(23) 


where A = h^j, — On the t-plane operator insertions 
are at U = 7r(2z + x)/n and = 7r(2z — x)ln, to = 
to = —TTx so (221 needs to be regulated. However, the 
regulators cancel in the expression for Jg'- 


Js = 


f nR 


sin(7ra;) 


-2(ti-1)(A-|-A) 


V 27r (cos(nt) - cos{-kx)) ) 

TD ’ / \ \ 2('n'—1)(A+A) 

R sin(7rj:) \ 

,27r (cos(t)-cos(7ra;))yj^^^ 




(24) 


Hence, in the limit n —>■ 1 there is no contribution to ([^ 
from the change of coordinates to the strip. A similar 
argument shows that the conformal factors do not con¬ 
tribute to the expression for the modular operator either. 
The same is true for using correlators on the full complex 


plane in (22). In other words, we have 


/sin(7r(z -|- 1/2)/R) \ 
\sin(7r(2; — l/2)/R)) 


/ sin(7r(z -|- l/2)/R)\ 
\sin(7r(z — l/2)/R)) 


( 20 ) 


uniformize the n-sheeted Riemann surface TZn to the com¬ 
plex plane and a strip of width 27r, respectively. For sim¬ 
plicity we focus on primary fields of dimension (h^,/i.^) 
which transform according to 


S{(j)\\tj;) = dn\og 




mH{^)-H{ao))\c/) 


dn^Og 


($$)(^f^)"-l 


(25) 


where the correlators can be either on the strip or the 
full complex plane, and we have suppressed the location 
of operator insertions. 




( 21 ) 


Then, 


i9z\ —Rn 


sin(7ra:) 


duj J 27ruj (cos(in log(uj)) — cos(7rx)) 
dz'\ inR sin(7ra:) 

Ift ^ 


( 22 ) 


27r (cos(nt) — cos(7ra:)) 

Under a change of coordinates ([^ transforms as 
Sit/H) = dn log [Js{{ti\, h^, h^)Js{{Q, h^, 

(oi”))5(4/(Fo)4/(to))r' 


+dn log 


($(t(,)4>(t„)oi”-i))s($(t')$(to))r' 


Analytic continuation 

We are interested in finding the analytic continuation 
of following finite sums 

n—1 n—2 n—1 

n n /n.m(x)"-'"-\ J] 

m—1 m—1 m—1 

where fn,m is given by 


fin,m) = - 


sin^ {nm/n) 


n—>-1 


sin(7r(m -|- a:)/n) sin(7r(m — x)/n) 


(26) 
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The answer is given by the identities below: 

n—1 


n 

Tl—l 
n—2 

n 

n—l 

n—1 


n—m— 1 


nsin(7rx/n) 

sin(7rx) 


n sin{Trx/n] 


sin(7ra;) 


n—2 


m—1 

n—l 




nsin(7rx/n) 
sin(7rx) 

sin(“+7)/3( 

7T7n/n) 


sin“^(7r(m + a;)/n) sin'*'^(7r(TO — x)/n) 

' nsin{TTx/n)'\^°'^^^^ 
sin(7ra:) 


( 27 ) 


To prove these identities, first bring the expression on 
the right to the left-hand side. Then, note that both 
the numerator and denominator always have poles of the 
same degree, and therefore has no poles at any value of x. 
Taking the limit x —>■ 0 we find that our four expressions 
are equal to one for all x, n, a, /3 and 7 . 



